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ABSTRACT 


A  general  computer  program  was  written  for  the 
Electrodata  204  "Datatron"  Digital  Computer  to 
calculate  acceleration  performance  and  related 
characteristics  of  tracked  vehicles. 

Vehicle  performance  is  determined  at  prescribed 
time  intervals  during  acceleration  from  stand¬ 
still  to  maximum  velocity.  Other  additional  data 
supplied  by  this  program  are  vehicle  speed,  dis¬ 
tance  traveled,  time  of  travel,  sprocket  torque, 
tractive  effort,  rolling  resistance,  drawbar  pull, 
and  acceleration. 
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PROJECT  TITLE:  DIGITAL  COMPUTER  PROGRAM  FOR  ACCELERATION 

PERFORMANCE,  TRACKED  VEHICLES 


INTRODUCTION: 

This  report  presents  a  high-speed  computer  program 
for  determining  maximum  acceleration  performance  of 
tracked  vehicles  on  level  or  grade  terrain. 

Data  resulting  from  this  program  are  vehicle  speed, 
distance  traveled,  time  of  travel,  sprocket  torque, 
tractive  effort,  rolling  resistance,  drawbar  pull, 
and  acceleration. 


OBJECT: 

Develop  a  digital  computer  program  for  rapid  and 
accurate  solution  of  vehicle  acceleration  perform¬ 
ance  . 


SUMMARY : 

1.  Study  of  vehicle  acceleration  is  a  specific 
application  of  Newton's  third  law  of  motion 
(F  =  MA) .  Synthesis  of  the  essential  force 
and  mass  factors  describes  acceleration  per¬ 
formance  of  track  type  vehicles  as  equal  to 
the  ratio  of  vehicle  drawbar  pull  and  equiv¬ 
alent  mass. 

F  =  MA 

A  =  F  z  Vehicle  Drawbar  Pull _ 

M  Vehicle  Mass  @  Acceleration 

-  Tractive  Effort  -  Resistance 
Equivalent  Mass 

I  [(tT  x  RFD  X  EFD)/PDsJ  -  twt  (Sine  0)  +  RrI 
(Wt/g)  x  (1.07  +  o0025  $  2) 

T^  =  Transmission  Torque 
Rpp  =  Final  Drive  Ratio 

Final  Drive  Efficiency 
Pitch  Diameter  of  Sprocket 
Vehicle  Weight 


EFD  = 
PDS  = 
Wt  = 


1 


0  -  Grade  Angle 

RR  -  Rolling  Resistance 

g  r  Gravity 

V  -  Gear  Reduction 

2.  Data  of  velocity,  time,  and  distance  are  auto¬ 
matically  plotted  on  graphs  of  time  vs.  velocity, 
time  vs.  distance,  and  distance  vs.  velocity. 

3.  Accuracy  of  the  program  was  established  by  com¬ 
paring  calculated  results  of  the  T113E1  with  actual 
field  tests.  Field  tests  and  calculated  accelera¬ 
tion  performance  are  illustrated  in  Figure  Al,  Ap¬ 
pendix  A. 

4.  Equipment  necessary  to  operate  the  program  is  the 
main  "Datatron"  computer  and  console,  a  magnetic 
tape  unit,  a  high-speed  punch  and/or  wide  carriage 
Flexowriter  with  format  control,  and  the  PACE  data¬ 
plotter. 

CONCLUSIONS  AND  RESULTS: 

1.  In  the  course  of  this  study,  it  was  necessary  to 
use  a  mathematical  model  representing  standard 
tracked  vehicle  design.  Unique  vehicle  designs 
may  require  program  modification  to  reflect  ade¬ 
quate  mathematical  precision. 

2.  Final  results  correlate  closely  to  actual  field 
tests  for  accurate  input  data.  The  program  offers 
an  opportunity  to  study  changes  in  over-all  vehicle 
performance  by  changing  characteristics  of  individ¬ 
ual  systems. 

3.  Accuracy  of  results  largely  depends  on  the  validity 
of  input  data.  Program  correlation  with  experi¬ 
mental  data  taken  in  the  field  indicates  accuracy 
within  limits  of  ±5$°  Program  solution  requires 
approximately  11  minutes  of  machine  time  and  18  min¬ 
utes  for  data  printout. 

BACKGROUND  INFORMATION: 

1.  Major  components  of  the  mathematical  model  presented 
in  Summary,  para.  1.,  other  than  those  of  basic 
vehicle  mechanics,  are  as  described  by  Jaroslav  J. 
Taborek  in  "Mechanics  of  Vehicles".1 


1.  Available  at  Reader  Service  Department,  MACHINE  DESIGN, 
Penton  Building,  Cleveland  3*  Ohio. 


2 


2.  Mr.  Taborek 1  s  inertia  resistance  formula,  =  1  + 
a  +  .0025  g  2jl”,  gives  r  as  an  equivalent  mass, 
a)  as  a  constant  reflecting  inertia  resistance  of 
the  track,  and  f  as  a  gear  reduction  between  engine 
and  sprocket.  The  equivalent  mass  is  a  variable 
that  takes  into  consideration  acceleration  of  rotat¬ 
ing  masses  within  the  vehicle. 


Page  39  of  "Mechanics  of  Vehicles" . 
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FLOW  CHART: 


The  flow  chart  (Figure  1)  is  a  graphic  explanation  of 
the  program.  The  program  sequence  and  operational  steps 
are  shown  together  with  the  method  of  using  input  data 
and  logic  of  decision  branching.  Figure  2  is  a  detailed 
representation  of  the  digital  computer  program  describ¬ 
ing  the  solution  of  vehicle  acceleration. 
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DIGITAL  COMPUTER  FLOW  CHART 
TRACKED  VEHICLE  ACCELERATION  PROGRAM 


•M 


-►C 


A  (START) 
FIGURE  2  (PART  2) 
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TEST  PROGRAM  INPUT  DATA: 


Test  program  input  data  are  illustrated  in  Appendix  B. 
Problem  solution  and  data  printout  require  the  use  of 
10  subroutines.  Subroutine  numbers  and  descriptions 
are  listed  in  the  same  order  as  they  appear  on  the  pro¬ 
gram  printout  (Appendix  C). 

4l0  Dataplotter  routine 

007  Automatic  high-speed  tape  feed 

400  Data  printout 

401  Linear  interpolation 

099  Error  routine 

204  Cosine  0 

203  Sine  0 

206  Arctan  0 

002  Alphanumeric  printout 

FPA.  Floating  point  arithmetic 


Compiler  code,  description  and  storage  locations  for 
program  variables  and  constants  are  as  follows  (these 
compiler  codes  are  shown  In  the  flow  chart  of  Figure  2) 


Compiler 

Notation 

Storage 

Location 

Description 

i 

1 

0001 

Variables  40-59  subscript 
for  grades  in  %  slope 

X2 

0002 

Variables  60-70  subscript 
for  shift  rpm 

i 

3 

0003 

Variables  80-89  subscript 
for  transmission  ratios 

i 

4 

0004 

Variables  90-99  subscript 
for  standard  transmission 
efficiencies 

i 

5  0005 
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Constant,  one  less  than  the 
number  of  gears,  for  stand¬ 
ard  transmission  only 


Compiler 

Notation 

Storage 

Location 

Description 

*6 

0006 

Variables  150-249  subscript 
for  Dataplotter  velocity 

i 

7 

0007 

Variables  250-349  subscript 
foo?  Dataplotter  time 

18 

0008 

Variables  350-449  subscript 
for  Dataplotter  distance 

19 

0009 

Subscript  for  Dataplotter 
index 

yo 

1 

0010 

1 

Alphanumeric  printout  of 
title  and  column  headings 

y99 

0109 

yioo 

0110 

Tan.  of  slope 

yioi 

0111 

Grade  resistance  *  lb  force 

y102 

0112 

Max.  tractive  effort 

yi03 

0113 

Transmission  rpm 

y104 

0114 

Velocity  mph 

yio5 

0115 

Engine  rpm 

y106 

0116 

Transmission  torque,  lb-ft 

yi07 

0117 

Engine  torque,  lb-ft 

yio8 

0118 

Rolling  resistance,  lb/ton 
and  lb  force 

yl09 

0119 

Sprocket  torque,  lb-ft 

yno 

0120 

Tractive  effort,  lb  force 

ym 

0121 

Equivalent  mass,  slugs 

yH2 

0122 

Acceleration,  ft/sec/sec 

10 

Compiler 

Notation 

Storage 

Location 

Description 

yH3 

0123 

A  t?  sec. 

yll4 

0124 

t,  sec. 

yll5 

0125 

A  d,  ft 

yil6 

0126 

d,  ft 

yH7 

0127 

Drawbar  pull,  lb  force 

yll8 

0128 

Velocity  printout  Index 

yl50 

| 

0160 

\ 

Storage  locations  for  Data¬ 
plotter  velocity 

y249 

0259 

y250 

0260 

1 

Storage  locations  for  Data¬ 
plotter  time 

y349 

0359 

y350 

i 

0360 

1 

Storage  locations  for  Data¬ 
plotter  distance 

y449 

? 

0459 

co 

0461 

Vehicle  weight,  lb 

C1 

0462 

Pitch  diameter  of  sprocket, 
in. 

c2 

0463 

Final  drive  ratio 

c3 

0464 

Final  drive  efficiency 

c4 

0465 

Coefficient  of  frictions,  % 

°5 

0466 

Max.  transmission  speed, 
rpm  (auto,  only) 

°6 

0467 

Max.  engine  speed,  rpm  (std 
trans.  only) 
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Compiler 

Storage 

Notation 

Location 

Description 

C7 

0468 

Engine  .shift  speed,  rpm 
(std.  trans.  only) 

°8 

0469 

Equiv.  mass  factors  for 
track  (1.07) 

c 

9 

0470 

A  V,  mph 

0 

1 — 1 

O 

0471 

Gear  shift  time,  sec  (std. 
trans.  only) 

c 

11 

0472, 

Trans,  decision  (1  for  auto 
0  for  std.) 

c 

12 

0473 

Printout  Interval,  mph 

C13 

0474 

Maximum  time,  sec. 

°14 

0475 

Maximum  velocity,  mph 

C15 

0476 

Maximum  distance,  ft 

ci 

Grade,  %  slope 

1 

0501-0520 

ci 

2 

0521-0540 

rpm  for  shift  to  next  gear 
(auto,  only) 

°80 

0541 

Transmission  ratio 

0i 

1 

89 

0550 

c 

90 

0551 

Std.  transmission  only, 
corres.  efficiency  for 
above  gear  ratios 

°99 

0560 

c 

100 

0561 

Std.  trans . --engine  speed. 

rpm 

Auto,  trans . --trans .  speed, 
rpm 


C119  0580 
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Compiler 

Notation 

Storage 

Location 

Description 

C120 

0581 

Std.  trans. --corres .  eng. 
torque,  lb -ft 

Auto,  trans . --corres .  trans 
torque,  lb -ft 

°139 

0600 

cl40 

l 

0601 

1 

Auto,  trans.  only--corres . 
eng.  speed 

C 

159 

! 

0620 

ci6o 

1 

0621 

1 

Velocity  points  on  rolling 
resistance  curve,  mph 

°179 

0640 

ci8o 

1 

o64i 

J 

Rolling  resistance  points, 
Ib/ton 

C199 

0660 
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FIGURE  Al.  DIGITAL  COMPUTER  PROGRAM  EVALUATION 

T113E1  CORRELATION 


APPENDIX  B 


COMPUTER  PROGRAM  INPUT  DATA  (T113E1) 


APPENDIX  B 


COMPUTER  PROGRAM  INPUT  DATA  (T113E1) 


Vehicle  weight 

Pitch  diameter  sprocket 


Final  drive  ratio 


Final  drive  efficiency 


23,200  lb 
19.61  in. 

4.31:1 
82.75  % 


Coefficient  of  friction 


Maximum  transmission  speed 
Equiv.  mass  factor  for  track 
Incremental  increase  in  velocity 
Trans,  decision  constant 


4,000  rpm 


•  5  mph 

1  (auto,  trans.) 


Printout  Interval 


1  mph 


Maximum  time 


Maximum  velocity 
Maximum  distance 


99  sec 

Constants  for 
99  mph  full  scale  de¬ 
flection  in  auto- 
9,900  ft  matic  plotting 
routine 


Shift  speeds: 


to  2nd 


750  rpm 


to  3rd 
to  4th 
to  5th 


1,070  rpm 
1,500  rpm 
2,060  rpm 


to  6th 


2,900  rpm 
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APPENDIX  B  (Cont'd.) 


Transmission  Ratios: 

1st 

5.296:1 

2nd 

3.810:1 

3rd 

2.690:1 

4th 

1.936:1 

5th 

1.390:1 

6th 

1.000:1 

Transmission  Speed  Corres.  Trans.  Torque 

Corres.  Engine 

(rpm) 

(lb-ft) 

(rpm) 

4,000 

194 

4,000 

2,901 

272 

2,900 

2,900 

301 

4,000 

2,06l 

388 

2,900 

2,060 

407 

4,000 

1,501 

534 

2,800 

1,500 

582 

4,000 

1,061 

752 

2,750 

1,060 

854 

4,000 

741 

1,067 

2,650 

740 

1,188 

4,000 

580 

1,494 

3,100 

431 

1,513 

2,300 

430 

1,513 

2,675 

300 

1,940 

2,580 

0 

3,763 

2,360 

Rolling  Resistance  Data: 

Rolling  Resistance 

Corres . 

Vehicle  Velocity 

(lb/ton) 

(mph) 

146.6 

40 

77.6 

5 

77.6 

18 
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APPENDIX  B 

COMPUTER  PROGRAM  INPUT  DATA  PRINTOUT  (T113E1) 


+0000000461 
+5523200000 
+5219610000 
+5143100000 
+5082750000 
+  5070000000 
+51+40000000 
+0000000000 
+0000000000 
+5110700000 
+5050000000 
+0000000000 
+51 10000000 
+5110000000 
+52990000 00 
+5299000000 
+5499000000 
+0000000501 
+0000000000 

+ oooooooooq _ _ 

+0000000521 

+5375000000 

+5410700000  Shift  rpm 
+5415000000 
+5420600000 
+5429000000 

+7070000000 _ 

+0000000541 

+5152960000 

+5138100000  Transmission 
+5126900000  ratios 

+5119360000 
+5113900000 

+5110000000 _ 

+0000000565 

+5440000000 

+5429010000  Transmission 

+5429000000  rpm 

+5420610000 

+5420600000 

+5415010000 

+5415000000 

+5410610000 

+5410600000 

+5374100000 


+0000000585 

+5319400000 

+5327200000  Transmission  Torque 
+5330100000  x  .97  ib.  ft. 
+5338800000 
+5340700000 

+5353400000 

+5358200000 

+5375200000 

+5385400000 

+541 0670000 
+5411880000 
+54t 4940000 
+5415130000 
+5415130000 
+54i 9400000 

+5437630000_ _ 

+0000000605 

+5440000000 

+5429000000  Engine  speed  rpm 

+5440000000 

+5429000000 

+5440000000 

+5428000000 

+5440000000 

+5427500000 

+5440000000 

+5426500000 

+5440000000 

+5431000000 

+5423000000 

+5426750000 

+5425800000 

+5423600000 _ 

+0000000654 

+0000000000 

+0000000000 

+0000000000 

+0000000000 

+5314660000  . 

+5277600000  R-H.  lb/ton 

+5277600000 _ 

+0000000634 

+0000000000 

+0000000000 
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+5374000000 

+5358oooooo 

+5343100000 

+5343000000 

+5330000000 

+0000000000 


+0000000000 
+0000000000 
+5240000000  MPH 

+5150000000 
+0000000000 
+0000080000 _ 
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APPENDIX  C 


COMPUTER  PROGRAM  PRINTOUT 


A  statement  and  program  printout  is  a  product  of  the 
Purdue  Compiler.  Each  compiler  statement  is  printed 
out,  and  the  computer  commands  necessary  to  perform 
the  statement  are  printed  below  it.  The  four  digit 
number  to  the  extreme  left  is  the  command  cell  loca¬ 
tion.  The  algebraic  sign  and  ten  digit  number  is  the 
computer  command. 


APPENDIX  C 


COMPUTER  PROGRAM  PRINTOUT  -  ACCELERATION 
PERFORMANCE  OF  TRACKED  VEHICLES 


REMARKS 


i  GCOO 
yO  00 10 

cO  0461 
sO  0662 
bO  0?4C 


410 

e 

03 

3520 

0557 

007 

e 

01 

3500 

0000 

400 

e 

02 

3460 

0519 

4oi 

e 

04 

3380 

0531 

099 

e 

01 

3360 

0460 

204 

e 

05 

3260 

0472 

203 

e 

05 

3160 

0467 

206 

e 

06 

3040 

0525 

002 

0 

01 

3020 

0536 

f.p.E 

i.e 

16 

3500 

0501 

2si1  = 

=40 

f 

0740  +  0000  64  7019 
0741  +  0000  02  0001 

1siO=561  f 


>  Computer  preparation 


0742  +  0000  64  7018 
074$  +  0000  02  0000 


0e"2e, 990010"  f 

0744  +  0000  64  7017 

0745  +  0100  21  5020 


J 


L Title  and  column  heading 
\  printout 


J 
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Osi6=150  f 


0746  +  0000  64  7016 
0747  +  0000  02  0006 

0si7=250  f 

0748  +  0000  64  7015 
0749  +  0000  02  0007 

0si8=350  f 

0750  +  0000  64  7014 
0751  +  0000  02  0008 
0752  +  0000  30  0760 
0753  +  0200  21  3060 
0754  +  0000  00  0350 
0755  +  0000  00  0250 
0756  +  0000  00  0150 
0757  +  0000  99  0010 
0758  +  0000  00  0561 

0759  +  0000  00  004c 

3s 12=60  f 

0760  +  0000  64  7019 
0761  +  0000  02  0002 

4si3=80  f 

0762  +  0000  64  7018 
0763  +  0000  02  0003 

5si4=90  f 

0764  +  0000  64  7017 
0765  +  0000  02  0004 

6syl04=0.f 

0766  +  0000  64  7016 
0767  +  0000  02  0114 

Osyl l8=0.f 

0768  +  0000  64  7016 
0769  +  0000  02  0128 


Computer  preparation 
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Osyl 1 4=0. f 

0770  +  0000  64  7016 
0771  +  0000  02  0124 


Osyl 16=0 of 

0772  +  0000  64  7016 
0773  +  0000  02  0126 

J 

7 syl  00="206e,  cil "  f  Find  o  from  tan  9 

0774  +  0000  72  0001 
0775  +  0000  30  0780 
0776  -  0000  00  0000 
0777  +  0000  00  0090 
0778  +  0000  00  0080 
0779  +  0000  00  0060 

0780  -  0000  64  0461 
0781  +  0100  21  3040 
0782  +  0000  02  0110 

8syl01=c0x("203e,yl00")  f  Grade  resistance 

0783  +  0000  64  0110 
0784  +  0100  21  3160 
0785  +  0461  21  3600 
0786  +  0000  02  0111 

9sy102=c0x("204e,y100")xc4  f  Max.  tractive  effort 

0787  +  0000  64  0465 
0788  +  0000  02  6000 
0789  +  0000  64  0110 
0790  +  0100  21  3260 
0791  +  6000  21  3600 

0792  +  0461  21  3600 
0793  +0000020112 

10syio3=(336.13445xylo4xc2)/c1  f  Transmission  rpm 

0794  +  0000  64  0462 
0795  +  0000  02  6000 
0796  +  0000  64  0463 
0797  +  0114  21  3600 
0798  +  0000  30  0800 
0799  +  0000  00  0060 
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0800  +  7019  21  3800 
0801  +  6000  21  3596 
0802  +  0000  02  0113 

1 1rgl2,rc1l/l  .f  Decision  for  auto,  or  std .  transmission 

0803  +  oooo  64  7018 
0804  +  0000  02  6000 
0805  +  0000  64  0472 
0806  +  6000  21  3584 
0807  +  0000  15  7011 
0808  +  0000  73  7017 
0809  +  0000  28  7011 
0810  +  0000  20  0674 

0g2lf  Go  statement 

0811  +  0000  2o  0683 


I2syl05=y103xci3  f  Engine  rpm 


0812  +  0000  72  0003 
0813  -  0000  64  0461 
081 4  +  0113  21  3600 
0815  +  0000  02  0115 
08 16  +  0000  30  0820 
08 17  -  0000  00  0000 
0818  +  5110  00  0000 
0819  +  5333  61  3445 


13rg16,xy105/=c7  f  Decision  for  transmission  shift 

0820  +  0000  64  0468 

0821  +  0000  02  6000 

0822  +  0000  64  0115 

0823  +  6000  21  3584 

0824  +  0000  15  7027 

0825  +  0000  73  7019 

0826  +  0000  28  7028 

0827  +  0000  2q  0678 


15si3=i3+1 


f  Index  std .  transmission 
ratio  to  next  gear 


0828  +  0000  30  0836 
0829  +  0000  74  0003 
0830  +  2000  29  3360 
0831  +  0000  02  0003 
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n  -l.i  i  Index  std.  transmission 
1  1  efficiency  to  next  value 

0632  +  0000  64  7018 
0833  +  0000  74  ooo4 
0634  +  2000  2 9  3360 
0835  +  0000  20  1308 

1 4rg1 9 , ri3/ (60-KL5 )  f  Test  for  last  transmission  ratio 

0836  +  0000  64  0005 
0837  +  0000  30  0840 
0838  +  0000  00  0001 
0839  -  0000  00  0000 

06 40  +  0000  74  7019 
0841  +  2000  29  3360 
0842  +  0000  02  6000 
0843  +  0000  64  0003 
0844  +  0000  75  6000 
0845  +  1000  29  3360 
0846  +  0000  15  7050 
0847  +  0000  73  7oi8 
0348  +  0000  28  7050 
0849  +  0000  2o  0681 

16rg69,ryi05/=c6  f  Test  for  max.  engine  speed 

0Q50  +  0000  64  0467 
0851  +  0000  02  6000 
0552  +  0000  64  0115 
0853  +  6000  21  3584 
0854  +  0000  15  1322 
0855  +  0000  73  7018 
0856  +  0000  28  1320 
0857  +  0000  30  0860 
0858  -  0000  00  0000 
0859  +  0000  00  0080 

0860  +  0000  20  0731 

0g42f  Go  statement 

0861  +  0000  20  0704 

I9sy105=y103xci3  f  New  engine  speed 

0862  +  0000  72  0003 
0863  -  0000  64  0461 
0864  +  0113  21  3600 
O865  +  0000  20  1305 
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1 7sy1 1 4=yl  1 4+clC  f  Time  adjust  while  shifting 


0866  +  0000  64  0471 
0867  +  0124  21  3580 
0868  +  0000  02  0124 


iSsyl 16=y1 16+(y104xe10) 


f  distance  adjust  while  shifting 


0869  +  0000  64  0471 
0870  +  0114  21  36OO 
0871  +  0126  21  3580 
0872  +  0000  02  0126 

69syi07="4oie,y105,0,1 

0873  +  0000  64  7019 
0874  +  0000  02  396O 
0875  +  0000  64  7018 
0876  +  0000  02  396I 
0877  +  0000  30  0880 
0878  -  0000  00  0000 
0879  +  0000  00  0001 
0880  +  0000  64  0115 
0881  +  0300  21  338O 
0882  +  0000  02  0117 

20sy 1 06=y1 07xci3xci4 


$  Find  engine  torque 
from  engine  rpm 


±‘  Transmission  torque 


0883  +  0000  72  0004 
0884  -  0000  64  0461 
0885  +  0000  72  0003 
0886  +  0461  21  3602 
0887  +  0117  21  3600 
0888  +  0000  02  0116 


0g24f  Go  statement 
0889  +  0000  20  0686 


21rg29,ryl03/=c5  f  Test  for  max.  transmission  rpm 


0890  +  0000  64  0466 
0891  +  0000  02  6000 
0892  +  0000  64  0113 
0893  +  6000  21  3584 
0894  +  0000  15  7097 
0895  +  0000  73  7019 
0896  +  0000  38  0900 
0897  +  0000  20  0691 
0898  +  0000  30  0900 
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0899  -  0000  00  0000 
0g42f  Go  statement 
0900  +  0000  20  0704 

29rg22,ry103/ci2  f  Decision  to  shift  gears 

0901  +  0000  72  0002 
0902  -  0000  64  0461 
0903  +  0000  02  6000 
0904  +  0000  64  0113 

0905  +  6000  21  3584 
0906  +  0000  15  7010 
0907  +  0000  73  7019 
0908  +  0000  28  7010 
0909  +  0000  20  0684 

30si2=i2+i  f  Index  shift  speeds 

0910  +  0000  64  701O 
0911  +  0000  74  0002 
0912  +  2000  25  3360 
0913  +  0000  02  0002 

31si3=i3+i  f  Index  transmission  ratio 

0914  +  0000  64  7016 
0915  +  0000  74  0003 
091 6  +  2000  29  3360 
0917  +  0000  30  0920 
0918  +  0000  00  0001 
0919  -  0000  00  0000 

0920  +  0000  02  0003 

22syio5="4oi e;yi03,o,2"  f  Eng.  speed 

0921  +  0000  64  7019 

0922  +  0000  02  5960 
0923  +  cooo  64  7018 
0924  +  0000  02  3961 
0925  +  0000  64  0113 
0926  +  0300  21  5380  . 

0927  +  0000  02  0115 

23syl06="4oie,y103,o,1"  f  Transmission  torque 

0928  +  0000  '64  7017 
0929  +  0000  02  396O 
0930  +  0000  64  7018 
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0931  +  0000  02  3961 

0932  +0000640113 
0933  +  0300  21  3380 
0934  +  0000  02  0116 

0935  +  0000  30  0940 

0936  +  2000  29  3505 

0937  +  0000  00  0001 
0938  -  0000  00  0000 
0939  +  0000  00  0002 

24sylo8="4oi e,ylo4,3j4"  f  Unit  rolling  resistance 

0940  +  0000  64  7019 
0941  +  0000  02  3960 
0942  +  0000  64  7018 
0943  +  0000  02  3961 
0944  +  0000  64  0114 
0945  +  0300  21  3380 
0946  +  0000  02  0118 

25syio8=(yioSxcO)/2ooo.  f  Rolling  resistance 

due  to  motion 

0947  +  0000  64  7017 
09^  +  0000  02  6000 
0949  +  0000  64  0461 
0950  +  0118  21  3600 
0951  +  6000  21  3596 

0952  +  0000  02  0118 

26syi o8=y1o8+y1 01  f  Total  rolling  resistance 

0953  +  0000  64  01 1 i 
0954  +  0118  21  3580 
0955  +  0000  02  0118 
0956  +  0000  30  0960 
0957  +  5420  00  0000 
0958  +  0000  00  0003 
0959  +  0000  00  0004 

27-syl 09=yl 06xc2xc3  f  Sprocket  torque 

0960  +  0000  64  0464 
0961  +  0463  21  3600 

0962  +  0116  21  3600 
0963  +  0000  02  0119 

28syil0=(y109x24.)/d  f  Tractive  effort 

0964  +  0000  64  0462 
0965  +  0000  02  6000 
0966  +  0000  64  7019 
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0967  +  0119  21  3600 

0966  +  6000  2i  3596 

096 9  +  0000  02  0120 

37rg32,iy110/=y102  f  Decision  for  printout 

0970  +  0000  64  0112 

0971  +  0000  02  6000 

0972  +  0000  64  0120 

0973  +  6000  2i  3584 

0974  +  0000  15  1325 

0975  +  0000  73  7018 

0976  +  0000  28  1523 

0977  +  0000  30  0980 

0978  -  0000  00  0000 

0979  +  5224  00  0000 

0980  +  0000  20  0694 

38sy1l0=y102  f  Index  for  next  printout 

0981  +  0000  64  0112 

0982  +  0000  02  0120 

32sy111=(c0x(c8+(.0025x(  c2xci3)p2)))/32.174o  f 


0983 

+ 

0000 

64 

7019 

0984 

+ 

0000 

02 

6000 

0985 

+ 

0000 

64 

7018 

0986 

+ 

0000 

02 

6003 

0987 

+ 

0000 

72 

0003 

0988 

0000 

64 

0461 

0969 

+ 

0463 

21 

3600 

0990 

+ 

0000 

72 

6005 

0991 

+ 

0000 

02 

6019 

0992 

+ 

0000 

64 

7017 

0993 

+ 

0000 

20 

7095 

0994 

+ 

60 19 

21 

3600 

0995 

+ 

0000 

22 

7094 

0996 

+ 

0000 

30 

1000 

0997 

+ 

5110 

00 

0000 

0998 

+ 

0000 

00 

0002 

0999 

+ 

5232 

17 

4ooo 

1000 

+ 

7019 

21 

3600 

1001 

+ 

0469 

21 

3580 

1002 

+ 

0461 

21 

3600 

1003 

+ 

6000 

21 

3596 

ioc-4 

+ 

0000 

02 

0121 

Equivalent 

mass 


29 


33sy1l7=yllo+(-ylo8)  f  Drawbar  pull 

1005  +  0000  64  0118 

1006  +  0000  02  6019 

1007  +  0000  65  6019 

1008  +0120  21  358O 

1009  +  0000  02  0127 

34sy  1 1 2=yi  1 7/yi  1 1  f  Acceleration 

1010  +  0000  64  0121 

1011  +  0000  02  6000 

1012  +  0000  64  0127 

1013  +  6000  21  3596 

1014  +  0000  02  0122 

35sy1  i 3=  ( 1 .4666667xc9)/y1 12  f 

1015  +  0000  64  0122 

1016  +  0000  02  6000 

1017  +  0000  64  0470 

1018  +  0000  30  1020 

1019  +  4825  00  0000 

1020  +  7019  2l  3600 

1021  +  6000  21  359c 

1022  +  0000  02  0123 

36sy115=(y104+(c9/2.))xy113x1. 4666667  ±'  A  d 

1023  +  0000  64  7019 

1024  +0123  21  3600 

1025  +  0000  02  6000 

1026  +  0000  64  7018 

1027  +  0000  02  6002 

1028  +  0000  64  0470 

1029  +  6002  21  3596 

1030  +0114  2i  3580 

1031  +  6000  21  3600 

1032  +  0000  02  0125 

39rg4o,rylo4/yl  18  f  Decision  for  printout 

1033  +  0000  64  0128 

1034  +  0000  02  6000 

1035  +  0000  64  0114 

1036  +  6000  2i  3584 

1037  +  0000  50  1040 

1038  +  5120  00  0000 

1039  +  511 4  66  6667 
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104c  +  oooo  15  7044 

1041  +  oooo  73  7019 

1042  +  oooo  26  7044 

1043  +  0000  20  0702 

4lsyil8=yll8+cl2  f  Index  printout  interval 

1044  +  oooo  64  0473 

1045  +  0128  2i  3580 

1046  +  0000  02  0126 


0g42f  Go  statement 


1047  +  oooo  20  0704 

40rg42,ryl  12/0.  f  Test  for  negative  acceleration 

1048  +  oooo  64  7019 

1049  +  0000  02  6000 

1050  +  0000  64  0122 

1051  +  6000  21  3584 

1052  +  0000  15  7056 

1053  *-  +  oooo  73  7019 

1054  +  OOOO  28  7056 

1055  yK  OOOO  20  0704 

Og46f  Go  statement 


1056  +  0000  20  0708 

Printout 

42e,'400e,y104,y1l4,y116,y103,y105,y109,y11C,y10S,yl17,yn2,ci 

5"  f 


1057  +  oooo  72  0003 

1058  +  0000  30  1060 

1059  -  OOOO  00  oooo 

1060  r  OOOO  64  0461 

1061  +  OOOO  02  396O 

1062  +  OOOO  64  0122 
1065  +  OOOO  02  396I 

1064  +  OOOO  64  0127 

1065  +  OOOO  02  3962 

1066  +  OOOO  64  0118 

1067  +  0000  02  3963 

1068  +  OOOO  64  0120 

1069  +  OOOO  02  3964 

1070  +0000640119 

1071  +  0000  02  3965 

1072  +  OOOO  64  0115 
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1073  +  0000  02  5566 

1074  +  oooo  64  0115 

1075  +  oooo  02  5967 

1076  +  0000  64  0126 

1077  +  OOOO  02  5968 

1078  +  OOOO  64  0124 

1079  +  oooo  30  1080 

1060  +  OOOO  02  5969 

1081  +  OOOO  64  0114 

1082  +1100  21  5460 

45syi6=y104  f  Velocity  storage 

1085  +  oooo  64  011  4 

1084  +  0000  72  0006 

1085  -  OOOO  02  0C10 

44syi7=yii4  f  Time  storage 

1086  +  OOOO  64  0124 

1067  +  0000  72  0007 

1088  -  OOOO  02  0010 

45syi8=yli6  f  Distance  storage 

1089  +  OOOO  64  01 2u 

1090  +  oooo  72  0008 

1091  -  OOOO  02  0010 

0^4 9f  Go  statement 

1092  +  OOOO  2c  0711 


46sy1 1 4-yl  1 4+yl  1 5  f  Time 

1093  +  OOOO  64  0125 

ioy4  +  0124  21  5580 

1095  +  oooo  02  0124 

47sy116=y1 16+yl  15  i*  Distance 

1096  +  OOOO  64  0125 

1097  +  0126  21  558C 

1098  +  OOOO  02  0126 

10 99  +  OOOO  30  1100 

48syl 04=y  1  04+c9  f  Velocity 

1100  +  OOOO  64  047c 

1101  +  0114  21  5560 
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1102  +  0000  02  0114 
OglOf  Go  statement 

1103  +  0000  20  0672 

49si6=i6+i  f  Index  velocity  storage 

1104  +  oooc  64  7019 

1105  +  0000  74  0006 

1106  +  2000  29  3360 

1107  +  0000  02  0006 

50si7=i7+i  f  Index  time  storage 

1108  +  0000  64  7019 

1109  +  0000  74  0007 

1110  +  2000  29  3360 
nil  +  0000  02  0007 

51si8=i8+i  f  Index  distance  storage 

1112  +  0000  64  7019 
1H3  +  0000  74  0008 

1114  +  2000  29  336O 

1115  +  0000  02  0008 


Osi9=i6  f  Index  storage 

1116  +  0000  64  0006 

1117  +  0000  02  0009 
1118  +  0000  30  1120 
1119  +  0000  00  0001 

52rg53,rc1l/l  .f  Transmission  decision 


1120  +  0000  64  7019 

1121  +  0000  02  6000 

1122  +  0000  64  0472 

1123  +  6000  21  3584 

1124  +  0000  15  7028 

1125  +  0000  73  7018 

1126  +  0000  28  7028 
1127  +  0000  20  0715 


54rg55,iy103/=c5 


fMax.  transmission 
printout  decision 


1128  +  0000  64  0466 

1129  +  0000  02  6000 

1130  +0000640113 

1131  +  6000  21  3584 


1132  +  0000  1 5  7035 

1133  +  0000  73  7018 

1134  +  oooo  28  7036 

1135  +  oooo  2o  0717 

Og56f  Go  statement 

1136  +  0000  2o  0718 

1137  +  0000  30  1 1 4o 

11 36  -  0000  00  0000 

1139  +5110  00  OOOO 

53^g55,ryl05/=c6  f  Max.  engine  printout  decision 

1140  +  0000  64  0467 

1141  +  0000  02  6000 

1142  +  0000  64  0115 

1143  +  6000  21  3584 

1144  +  0000  15  7047 

1145  +  OOOO  73  7019 

1146  +  OOOO  28  7048 

1147  +  oooo  2o  0717 

Og56f  Go  statement 

1148  +  oooo  20  0718 

55rg56,ryll2^o.  f  Test  for  negative  acceleration 

1149  +  0000  64  7019 

1150  +  0000  02  6000 

1151  +  0000  64  0122 

1152  +  6000  21  3584 

1153  +  oooo  15  7057 

1154  +  oooo  73  7019 

1155  +  oooo  28  7057 

1156  +  0000  20  0718 

0g46f  Go  statement 

1157  +  OOOO  20  0708 

1158  +  OOOO  30  1160 

1159  -  OOOO  00  OOOO 

56e"7e,500"  f  Automatic  tape  feed 

1160  +  OOOO  64  7019 

1161  +  0100  21  3500 

57si6=i6+(-l)  f  Index  plotter  iteration 
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UN 


1162  +  0000  64  7018 
1165  +  0000  02  6019 

1164  +  0000  65  6019 

1165  +  0000  74  0006 

1166  +  2000  29  5360 
1167  +  0000  02  0006 

8si7  =i7  +(-l)  f  Index  plotter  iteration 

1168  +  0000  64  7018 

1169  +  0000  02  6019 

1170  +  0000  65  6019 

1171  +  0000  74  0007 

1172  +  2000  29  5360 

1173  +  0000  02  0007 

60e"4l0e,0,yi7,yi6,cl5,cl4"  f  Plotter  printout 

1174  +  0000  64  0475 

1175  +  0000  02  596O 

1176  +  0000  64  0474 

1177  +  0000  50  1180 

1178  +  0000  00  0001 

1179  +  0000  00  0500 

1180  +  0000  02  5961 

1181  +  0000  72  0006 
1182  -  0000  64  0010 
1185  +  0000  02  5962 

1184  +  0000  72  0007 

1185  -  0000  64  0010 
1166  +  0000  02  3955 

1187  +  oooc  64  7019 

1188  +  0500  21  5520 

65rg71  ,rio/l51f  Test  further  iterations 

1189  +  0000  64  7018 

1190  +  0000  02  6000 

1191  +  0000  64  0006 

1192  +  0000  75  6000 

1193  +1000  29  356O 

1194  +  0000  15  1526 

1195  +  0000  73  7019 

1196  +  0000  28  1526 

1197  +  0000  50  1200 

1198  +0000000151 

1199  -  0000  00  0000 

1200  +  OOOC  20  0755 
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Og57f  Go  statement 

1201  +  0000  20  0719 

71e"7e,500"  f  Automatic  tape  feed 

1202  +  0000  64  7019 

1203  +  0100  21  3500 

0si8=350  f 

1204  +  0000  64  7018 

120$  +  0000  02  0008 

74e"4l0e,0,yi8,yi6,c!5Jcl2fn  f  Plotter  printout 

1206  +  0000  64  047$ 

1207  +  0000  02  3960 

1208  +  0000  64  0476 

1209  +  0000  02  3961 
1210  +  0000  72  0006 

1211  -  0000  64  0010 

1212  +  0000  02  3962 
1213  +  0000  72  0008 

121 4  -  0000  64  0010 

1215  +  0000  02  3963 

1216  +  0000  30  1220 

1217  +  0000  30  1200 

121 8  +  0000  00  0350 
1219  +  0000  00  0500 

1220  +  0000  64  7019 

1221  +  0500  21  3520 

72si6=i6+l  f  Index  plotter  iteration 

1222  +  0000  64  7018 

1223  +  0000  74  0006 

1224  +  2000  29  3360 

1225  +  0000  02  0006 

73si8=i8+l  t  Index  plotter  iteration 

1226  +  0000  64  7018 

1227  +  0000  74  0008 

1228  +  2000  29  336O 

1229  +  0000  02  0008 

75rg65,ri6z(i9)+1  f  Test  for  further  iterations 
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1250  +  0000  64  7018 

1231  +  0000  02  o000 

1232  +  0000  64  0009 

1233  +  0000  74  6000 

1234  +  5000  29  336O 

1235  +  0000  02  6000 

1236  +  0000  64  0006 

1237  +  0000  30  1240 

1238  +  0000  00  0001 

1230  _  0000  00  0000 

1240  +  0000  75  6000 

1241  +  1000  29  3360 

1242  +  0000  04  7044 

1243  +  0000  20  0727 

0g7’4f  Go  statement 

1244  +  0000  20  0736 

65e"7e, 500"  f  Automatic  tape  feed 

1245  +  0000  64  7019 

1246  +  0100  21  3500 

0si8=350  f  Computer  preparation 

1247  +  0000  64  7018 

1248  +  0000  02  0008 


61  e"4l0e,0,yi7,yi8,c13,c15"  f'  Plotter  printout 

1249  +  0000  64  0476 

1250  +  0000  02  396O 

1251  +  0000  64  0474 

1252  +  0000  02  396I 

1253  +  0000  72  0008 

1254  -  0000  64  0010 

1255  +  0000  02  3962 

1256  +  0000  72  0007 

1257  +  0000  30  1260 

1258  +  0000  00  0350 

1259  +  0000  00  0500 

1260  -  0000  64  0010 

1261  +  0000  02  3963 

1262  +  0000  64  7019 

1263  +  0500  21  3520 

0si8=i8+1f  Index  plotter  iteration 
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1264  +  0000  64  7018 

1265  +  0000  74  0008 

1266  +  2000  29  3560 
1267  +  0000  02  0008 

Osi7=i7+1f  Index  plotter  iteration 

1268  +  0000  64  7018 

1269  +  0000  74  0007 

1270  +  2000  2q  3360 

1271  +  0000  02  0007 

0rg64,ri7zi9+101  f  Test  for  further  iteration 

1272  +  0000  64  7017 

1273  +0000  74' 0009 

1274  +  2000  29  336O 

1275  +  0000  02  6000 

1276  +  0000  30  1280 

1277  +  0000  00  0101 

1278  +  0000  00  0001 

1279  _  0000  00  0000 

1280  +  0000  64  0007 

1281  +  0000  75  6000 

1282  +  1000  29  3360 

1283  +  0000  o4  7085 

1284  +  0000  20  0726 

Ogolf  Go  statement 

1285  +  0000  20  0723 

64e"7e,500"  f  Automatic  tape  feed 

1286  +  0000  64  7019 

1287  +  0100  21  5500 

66sil=il+l  f  Index  slope  subscripts 

1288  +  0000  64  7018 

1289  +  0000  74  0001 

1290  +  2000  29  3360 

1291  +  0000  02  0001 

67rg68,rci1z0.f  Test  for  last  run 

1292  +  0000  64  7017 

1293  +  0000  02  6000 

1294  +  0000  72  0001 
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0000 

30 

0866 

1310 

+ 

0000 

00 

0000 

1311 

+ 
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